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Abstract 

Data collected with the BABAR detector at the PEP-II asymmetric B Factory at SLAC are used 
to measure the lifetimes of the B° and B + mesons. The data sample consists of 8.3 fb _1 collected 
near the T(4S) resonance. B° and B + mesons are fully reconstructed in several exclusive hadronic 
decay modes to charm and charmonium final states. The B lifetimes are determined from the flight 
length difference between the two B mesons. The preliminary measurements of the lifetimes are 

t b o = 1.506 ± 0.052 (stat) ± 0.029 (syst) ps, 
t b+ = 1.602 ± 0.049 (stat) ± 0.035 (syst) ps, 

and of their ratio is 

t b +/t b o = 1.065 ± 0.044 (stat) ± 0.021 (syst). 
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1 Introduction 



This paper presents preliminary measurements of the and B°/B° lifetimes and their ratio 
with data collected by the BABAR detector jjj at the PEP-II asymmetric-energy B Factory. At 
the T(4S), B mesons are produced in B + B~ or B°B° pairs. In PEP-II, the center-of-mass frame 
is boosted roughly along the z axisQ (/?7 ~ 0.56). In this analysis, one of the B mesons is fully 
reconstructed in a variety of clean hadronic decay modes to charm and charmonium final states. 
An inclusive technique is performed to reconstruct the decay vertex of the opposite B meson, and 
the lifetime is determined from the distance along the z-axis between the decay vertices of the two 
B mesons (Az = z rcc — z opp ). The true value of Az is distributed exponentially with an average of 
(\Az\) = (p-f) z B cT B ~ (Pr (4S) /m r(4 s))cr B . 

2 Event sample and vertex reconstruction 

The data used in this analysis were collected by the BABAR detector at the PEP-II storage ring 
in the period from January to June 2000. The total luminosity is 7.4 fb -1 collected at the T(4S) 
resonance and 0.9 fb -1 40 MeV below the resonance. The number of produced BB pairs is estimated 
to be 8.4 x 10 6 . 

B° and B + mesons are reconstructed in the following modes (and their charge conjugates): 
B ° _» DW-7T+ D^-p + , D^'af, J/^K*° and B + -» D ( * )0 ir+, J/i;K+, ip(2S)K+. All final 
state particles are reconstructed. The signal region for each decay mode in the selected sample is 
defined by the three standard deviation bands in the two-dimensional distribution of the kinematic 
variables, AE and mES0- The resolution on m-ES is about 3MeV/c 2 , and that on AE varies from 
mode to mode between 12 and 40 MeV. The purities of final samples are approximately 90% 
(Fig.l). 

A geometric and kinematic fit of the fully reconstructed B Tec is performed and is required to 
converge. The vertex resolution ranges from 45 to 65 /im, depending on the mode. The vertex 
of the opposite B is determined using all the tracks that are not associated with the B Tec , along 
with the -B pp "pseudotrack" , which is estimated from the momentum difference between T(4S) 
and -B rec , and the beamspot. The fitting procedure is repeated, after removing tracks and V° 
decays that result in poor fits, until no track contributes more than a certain value to \ 2 - This 
algorithm gives one standard deviation errors of 115 //m and biases around 25 //m with Monte 
Carlo simulation. The Monte Carlo simulation is also used to study the Az resolution function of 
the residual 5(Az) = (^z)mcasured — (<^<z)generated and the pull 5(Az)/a(Az). Two-Gaussian fits 
successfully describe the Monte Carlo residual and pull distributions. They give a one standard 
deviation width of 130 //m and a bias of 24 /im for the residual. Because it is dominated by the 
-Bopp vertex, the resolution function shape is essentially the same for all modes we have considered. 

3 Lifetime fits and systematic uncertainties 

The lifetime is extracted from the Az distribution of the selected events with an unbinned maximum 
likelihood fit. The measurements performed on each event i are represented by three input numbers: 

a The axis of the PEP-II beams is tilted by 20 mrad with respect to the BABAR z axis. 

h AE = E* ec — El, m-ES = \J E"^ — p*£ c where E* ec and p* oc are the measured B candidate energy and momentum, 
and El is the beam energy, all denned in the center-of-mass frame. 
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(Az)i, its error cjj, and the probability for an event i to come from signal based on the tbes spectrum. 
All events that satisfy 5.2 < tties < 5.3GeV/c 2 and have AE in the signal region are input to the 
fit. The fit determines the B lifetime, the proportion of outliers, and two sets of parameters that 
describe the resolution function and the Az shape of the background. The signal distribution is 
represented by a convolution of the theoretical Az distribution (two exponential wings) and the 
pull representation of the Az resolution function. A function of the form G <g> (1 + E) is used[]. 
The background Az distribution is described by the sum of a single Gaussian and two independent 
exponential tails, one for positive Az and one for negative Az. The outliers are described by a 
wide symmetric Gaussian with a fixed width of 2500 /jm. 

We fit the Az distributions of B° and B + samples simultaneously with a single set of parameters 
for the resolution function but different sets of parameters (one per charge) to describe the lifetime, 
the background and the outliers (Fig.l). The lifetime ratio is determined with a similar combined 
fit with the B° lifetime and the ratio being free parameters. The statistical errors are about 
3.5%. About 2% of the statistical error is due to the high correlation between lifetime and the Az 
resolution function parameters. 

The dominant systematic uncertainties are the limited Monte Carlo statistics for checking Az 
distortion due to event selection, the modeling of Az outliers, the length scale measured along 
z, and the background modeling. Each contributes about 1% to the systematic uncertainty. All 
other systematic uncertainties we have studied, including the choice of parameterization for the 
resolution function and errors on the boost measurements, are less than 0.5%. The approximation 
of using the T(4S) boost for both B mesons in an event results in a 0.4% shift. The central values 
of the lifetimes are corrected for this shift. 




Beam-Energy Substituted Mass (GeV/c) Deca y time difference A t (ps) 



Figure 1: Left: toes distribution for all the hadronic modes for B , fit to the sum of a Gaussian 
distribution and the ARGUS background parameterization]^]. The total yields in all B° and B + 
modes are 2210 ± 58 and 2261 =L 53 (not shown), respectively. Right: At distribution for B° /B°. 
The result of the lifetime fit is superimposed. The background is shown by the hatched distribution. 



c The sum of an unbiased Gaussian and the convolution of the same Gaussian with a decaying exponential. This 
resolution function leads to the smallest overall error with present statistics. 
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The preliminary results for the B meson lifetimes are 

t b o = 1.506 ± 0.052 ± 0.029 ps, t b+ = 1.602 ± 0.049 ± 0.035 ps, 
and for their ratio is 

t b+ / Tb0 = 1.065 ± 0.044 ± 0.021. 

The first error is statistical and the second is systematic. The results are consistent with previous 
B lifetime measurements || and competitive with the most precise ones. 
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